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We show that the Amontons' law of friction breaks down in the elastic object which has the 
friction force obeying the Amontons' law locally. By the finite element method, it is shown that 
the bulk static friction coefficient /xm decreases with increasing pressure and system length. It also 
depends on the apparent contact area. These behaviors indicate the breakdown of the Amontons' 
law and result from the local slip before the onset of bulk sliding. These results are consistent with 
experiments. An effective one-dimensional model well explains these behaviors analytically. 

PACS numbers: 46.55.+d, 46.50.+a, 62.20.Qp, 81.40.Pq 



In the 15th century da Vinci discovered that the re- 
sistance against the sliding motion of the solid object 
on substrate does not depend on the apparent contact 
area and is doubled if the load of the object is dou- 
bled [U [2]. The discovery was rediscovered by Amon- 
tons about 200 years later and is now summarized by, 
the so-called, Amontons' law of friction: the friction force 
does not depend on the apparent contact area and is pro- 
portional to the loading force. The ratio of the friction 
force to the loading force is called the friction coefficient, 
which does not depend on the apparent contact area and 
the loading force. It is well known that the Amontons' 
law holds well for various systems. The law is usually 
explained as resulting from increase of the area of actual 
contact points, in which two surfaces actually contact 
each other via intermolecular forces, in proportion to the 
loading force [TJ |2] . 

It is also well known, however, that the Amontons' 
law breaks down in some systems, e.g., the systems with 
large adhesion [2j. Here, we show the breakdown of 
the Amontons' law in an elastic object by both of nu- 
merical and analytic methods. The object has the fric- 
tion force obeying the Amontons' law locally and is under 
uniform pressure Pext and pushed by the external shear 
force Ft from the trailing edge as shown in Fig. [ija). 
Numerical analysis based on a 2D finite element method 
(FEM) shows that the bulk static friction coefficient /iM, 
which corresponds to the onset of sliding of the whole 
object, decreases with increasing Pext- That indicates 
the breakdown of the Amontons' law and is consistent 
with experiments [3l [4]. The breakdown is closely re- 
lated to the local slip also observed in experiments [5]-[8] 
and theoretical and numerical works j8HT3]. The shear 
force Ft well below the bulk maximum static friction 
force induces a precursor, that is the local slip from the 
trailing edge of the object. The slip front proceeds slowly 
with slowly increasing Ft. When the traveling length of 
the front reaches a critical length the slow motion 
becomes unstable and changes to rapid motion. Then, 
the slip front quickly arrives at the trailing edge and the 
whole object slides. The magnitude of jiyi depends on 
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FIG. 1. (color online), (a) The elastic object is under uniform 
pressure Pext = Fn / {LW) and pushed at height /i by a rigid 
rod with the shear force Ft- (b) The ratio of the shear force to 
the loading force Ft/Fn vs the displacement of the pushing 
rod U = Vi for L = 1.0, d = 0.5, Pext = 0.0004 and the 
parameter set (A), (c) The local slip region is shown by the 
gray area (the red area for online version) in the U — x plane 
for the same parameters with those in (b). The horizontal 
line shows the critical traveling length of the slip front 

4 as m - {l^s - m){^c/L), where /is,/^K, and L 

are the local static and kinetic friction coefficients and 
the system length, respectively and /is > Mk- The crit- 
ical length £c decreases with increasing Pext- Then, fiu 
decreases with increasing Pext- These behaviors are well 
explained analytically based on a ID effective model. 

We first analyze the rectangular elastic object on a 
rigid substrate by FEM, which has length F, width W 
and height H along x, z-axes, respectively and mass 
density p. The plain stress condition is assumed along the 
axis. We apply the uniform pressure Pext = Fn/{LW) 
at the top surface, where Fn is the loading force, and 
the shear force Ft from the trailing edge at height h by 
a rigid rod. Hereafter, we set 1^ = 1 and h = H/2. 
The equation of motion of the displacement vector of 
the elastic object u{r^t) is expressed as p^u{r^t) = 
V • a [14 . Here, a = cr^®^^ + cr(dis) jg stress tensor 
composed from the elastic part cr^®^^ and the dissipation 
part t7*^^^^\ of which (a,/3) component is expressed as 

^S^^ = Vi^a(3^r]2{^xx^^zz)Sa(3' Here, ea(3 is the compo- 



nent of the strain rate tensor and rji , r]2 are the viscosity 
coefficients. The boundary condition at the top surface 
is given by azz = —Pe^t and axz = 0. The stress is free 
at X = and L except at the rigid rod, which has width 
of O.IH and moves with constant velocity V. The local 
friction force obeys the Amontons' law and gives the fric- 
tional stress, cr^^^^ = iJ. (ux) sign{ux)\azz\0{—crzz), at the 
interface, where 0{x) is the Heaviside function JTE\. The 
velocity v dependent friction coefficient jii{v) takes value 
from to /is at = 0, decreases linearly from jus to juk 
for < I'l'l < and equals to juk for Vc < \v\. 

We introduce the normalized quantities in which the 
length, mass and time are normalized by Lq = r]i/^/pE^ 
rriQ = r]i/\/ pE^ and to =^i/E^ respectively, where E is 
the Young's modulus. We express normalized quantities 
with tilde. The model has parameters L = L/Lo^d = 
H/L^Pq^i = Pq^i/E^V = Vto/Lo, which characterize 
the sample size and the external condition, /is7MK,^2 = 
^2/'^i,'^c7 and the Poisson's ratio which character- 
ize the material properties. We employ two parameter 
sets: (A) simulates a model viscoelastic material, which 
has u = 0.34, /is = 0.38, /ix = 0.1, Vc = 3.4 x lO""^, 
and 7)2 = 1; (B) simulates polymethylmethacrylate 
(PMMA), which has u = 0.4, /is = 1.2, /ix 0.2, Vc = 
3.9 X 10"'^, and 7)2 = 1 [16^- In the FEM calculation, 
the object is divided into equal-size cells with rectangu- 
lar shape. The number of cells is 40 x 40 or 80 x 80. The 
convergence of the results with respect to the number is 
verified. The following results are obtained for the limit 
of vanishing driving velocity. 

In Fig.jljb), we show the ratio of the shear force to the 
loading force Ft/Fn as a function of the displacement of 
the pushing rod U = Vt. At the beginning, we apply 
the uniform pressure at the top surface and make the 
elastic object relax. Then, we start pushing the object 
at t = 0. Periodic stick-slip motions are observed except 
the 1-st one, which is strongly infiuenced by the initial 
condition. Hereafter, we focus on the periodic regime. 
We define the bulk static friction coefficient /iM by the 
peak value of Ft/F^. The pressure Pgxt dependence of 
/iM is shown in Figs. [2|a,b) for the parameter set (A) 
and (B), respectively. It is to be noted that for PMMA, 
simulated by the parameter set (B), Pext = 1.0 x 10~^ 
and L = 0.02 correspond to Pext = 3.0 x 10^ [Pa] and 
L = 0.96[m], respectively [17]. Clear decreases of /iM 
with increasing Pext and L are observed, which indicate 
the breakdowns of the Amontons' law and the extensive 
property of the friction force. Such pressure dependence 
is consistent with experiments of PMMA |4]. It is 
well known that adhesion makes the friction coefficient 
increase with decreasing pressure [1] [2] , but it does not 
exist in the present system. We show Pn dependence of 
/iM for the case laying the longer edge of a rectangular 
object on the substrate and that laying the shorter one in 
Fig. |2|c). The friction coefficient decreases with increas- 
ing Pn and depends on the apparent contact area, which 




FIG. 2. (color online). The pressure Pext dependence of /iM 
for the parameter set (A) (a) and (B) (b). The lines with 
symbols indicate the results of the FEM calculation. Thin 
lines indicate the analytic results based on the effective ID 
model, where we set = 0.2. The lines with the same type 
correspond to the same value of L. (c) The loading force 
Pn dependence of /xm for the case laying the longer edge of 
a rectangular object on the substrate (L = 1.0, P = 0.5) 
and that laying the shorter one (L = 0.5, P = 1.0) for the 
parameter set fA). 
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FIG. 3. (color online). The bulk friction coefficient /iM as 
a function of ic/L for various values of Pext and L and for 
the parameter set (A) and (B). The lines show the theoretical 
result, Eq. (jsj. The inset shows the pressure Pext dependence 
of ic/L for d = 0.5 and the parameter set (A). 

also indicates the breakdown of the Amontons' law. 

The change of /iM shown in Fig. [2] results from the lo- 
cal slip before the onset of bulk sliding also observed in 
experiments of PMMA ^-^Sj. In Fig. [ijc), we show the 
local slip region, where the velocity at the interface has 
a finite value, in the U — x plane. The shear force Pt 
well below the bulk maximum static friction force induces 
slow local slip from the trailing edge. The slip front pro- 
ceeds slowly and the traveling length of the front i grows 
with U. When i reaches the critical length ic, the front 
starts rapid motion, which accompanies rapid local slip. 
Then, the front enters the leading edge quickly, the whole 
object slides, and Pt drops sharply. These behaviors are 
consistent with experiments [6]. Further increase of U 
repeats the similar behavior. Figure [3] shows that /iM is 
scaled by ic/L for each parameter set. The scaling with 
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FIG. 4. (color online) The normalized pressure p(x)/Pext (a) 
and the ratio of the shear stress to the pressure f{x)/p{x) (b) 
at the interface for the magnitudes of U indicated by arrows 
in Figs. [T]^b,c). The arrows in (b) indicate the position of the 
slip front. The horizontal lines indicate /is, Mm, and /xk- The 
parameters are same with those in Figs. [TJb,c). 

the pressure Pext dependence of ic/L shown in the inset 
of Fig. [3] gives the Pext dependence of jj^m- 

We show the pressure at the interface p{x) = —azz\z=o 
in Fig. |4[a ) for the magnitudes of U indicated by arrows 
in Figs7TTla,b). The applied pressure at the top surface 
is uniform, but p{x) depends strongly on x and increases 
with X. Similar pressure distribution is observed in ex- 
periments m |7] and results from the torque induced by 
the shear force Ft El- The low pressure at the trail- 
ing edge makes the local maximum static friction force 
small. As a result, the local slip starts from there. Fig- 
uie^h) shows the ratio of the shear stress at the interface 
f{x) = axzlz^o to the pressure, f{x)/p{x). Just after the 
sliding of the whole object, which accompanies rapid slip 
in the entire interface, the ratio is approximately equals 
to /iK (dashed line). The reason is considered as follows. 
During the acceleration of the rapid slip of any portion 
in the interface, the local friction coefficient there de- 
creases to /iK- At the instance of vanishing acceleration, 
the local shear stress at the interface f{x) is equal to the 
local frictional stress iiy^pix). Then, the portion deceler- 
ate and the local friction coefficient rapidly approaches 
to /is- When the portion stops slip, however, f{x) keeps 
the magnitude /j.kp{x) due to finite relaxation time. Fi- 
nally, the local frictional stress after the rapid sliding of 
the whole object relaxes to /j^kp{x). By contrast, in the 
region where the slow local slip has occurred, the fig- 
ure shows that the ratio f{x)/p{x) equals to /is because 
f{x) is given by the local frictional stress for vanishing 
velocity /isp(^) after the slow slip [18]. Figure |4]^b) also 
shows that the local magnitude of f{x)/p{x) can be larger 
than the bulk frictional coefficient /iM without precipitat- 
ing the local slip as observed in experiments jH |7] . The 
bulk maximum static friction coefficient jum is then de- 
termined at the onset of rapid motion of the local slip 
front, which is also shown analytically below. 

In order to analyze the transition from the slow to 
rapid motions of the local slip front at we employ an 
effective ID model, in which the degrees of freedom of 
the elastic object along the z-axis are neglected. The 



equation of motion of the model is expressed as, 
d'^u{x,i) ddxx{x,t) dxz{x,t) - d^xT'' 



dp 



dx 



OH 
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Here, u{x^t) is the displacement along the x-axis at the 
interface, and OH is the characteristic length of the vari- 
ation of the {xz) component of the stress. Here, we em- 
ploy ^ as a fitting parameter. The shear stress is given by 
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-{V-f~u), where 

rjt = 1 + 7)2 - The frictional stress axl^"^^ is given as before 
and El = !/{(! + - u)}, E2 = 1/{2(1 + u)}. We as- 
sume that the pressure at the interface is given by p(x) = 
2Pext^/^, which simulates the result shown in Fig. [4[a). 
The boundary condition is given by -^u{x^i) = at 
X = and L. Hereafter, we set the origin of U to the 
position of the pushing rod just after the sliding of the 
whole object. The adiabatic solution of Eq. ([T]) for [/ = 
is obtained as 



^ 2/iK/iPext r / u 

uq[x) = — — < ^(cosh 

^ ^ sinh 
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where = E1/E2 and = OhH. Here, we set the local 
friction coefficient to /ix as discussed before. For /7 > 0, 
the local slip front proceeds to £. The adiabatic solution 
for U > is expressed as 



Ua{x) 



ui{x), 
uo{x), 



0<x <£, 

i<x. 
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The value of Ua,{x) and its derivative must be continuous 
at X = ^. In the region I < x^ the increase of the shear 
stress with increasing U is cancelled by the increase of the 
local frictional stress, which takes value between fiKp{x) 
and /isp(x). The solution is then given by uo{x). The 
solution in the region < x < where the local friction 
coefficient is equal to /is, is obtained as. 



X — £ 

ui(x) = uo(x) -\-U(l — cosh ) 

, 2(/is - m)hPe^t r . . £-x 
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The boundary condition at x = gives £ 
E2LU /{{fis - I^K)hPe^t} for £/L <C 1. Hence, £ grows 
adiabatically with adiabatic increase of U as shown in 
the FEM calculation Fig. [ij^c). Similar result is obtained 
based on a ID spring-block model [T3^. When £ reaches 
the the critical length £c^ the adiabatic motion becomes 
unstable and changes to rapid motion. Then, the whole 
object slides. Employing the adiabatic equation, substi- 
tuting the relation between U and ^, and setting £ = £c 

in the expression of the shear force Pp = Jq dxdxz {x) , 
we get 



/^M = m + (/is - /iK)y- 
ij 



(3) 
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for Ic/L <C 1. This relation agrees well with the FEM 
calculation as shown in Fig. [s] even for Zc/Z < 1. 

The critical length 1^ is given by the linear stability 
analysis of Eq. ([T]) in the limit of vanishing driving ve- 
locity V. Setting u{x^i) = Ua{x) + Su{x^i) and lineariz- 
ing Eq. ([T]) with respect to 5u{x^i), we obtain its time 
evolution equation. The boundary condition is given 
by -^6u{x,i)\x=o = and Su{x > = 0. We ex- 
press du{x^i) = ^^^Q ii^ exp (cj^t) cos(/c^x/£), where 
km = (2m + l)7r/2. The time evolution equation gives 
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Here, off-diagonal terms are neglected, which are higher 
order of ^. The eigenvalue ujn gives the instability con- 
dition of Ug,. Usually, the instability results from the 
appearance of positive real part of any eigenvalue. In the 
present case, however, the backward motion of the fluc- 
tuation 5u is inhibited because it relaxes the local shear 
stress and the frictional stress stops it. That means the 
oscillating fluctuation never grows. The appearance of 
positive eigenvalue results in the instability. For the case 
of small value of Pext or ^, the viscosity in the 2nd term in 
Eq. Q stabilizes the adiabatic motion, but in the oppo- 
site case the velocity weakening friction force expressed 
by the last term leads to the instability. The instability 
occurs in the n = mode at first and gives the relation 
between the critical length 1^ and Pext as 
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This equation gives ic/L decreasing with increasing Pext 
and L: for 4/Z <C 1, Ic/L oc (LPext)"^/^- That is con- 
sistent with the FEM calculation shown in the inset of 
Fig. [sj The function 4 (Pext , L) is given by Eq. ([5| implic- 
itly. By putting into Eq. ([3|, we obtain the pressure 
and system length dependences of the friction coefficient 
fiM- The results are shown in Figs. [2|a,b). The analyt- 
ical results agree with the FEM calculation semiquanti- 
tatively for parameter set (A) and qualitatively for (B). 

In conclusion, we show that the bulk static friction co- 
efficient /iM of the elastic object decreases with increasing 
pressure and system length, which is consistent with ex- 
periments O [4| . The magnitude of fiu also depends on 
the apparent contact area. These behaviors have strong 
relation with the transition from slow to rapid motions 
of the local slip front. For large or small enough pres- 
sure and/or system length, where the ratio of the critical 



length of the transition to the system length is much less 
than or approximately equals to unity, fiu becomes con- 
stant and the Amonton's law holds. The present result 
may give a new method to control friction. 
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